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'The 3-phase version of the hybrid chiral bag model, containing the phase of asymptotic freedom, 
(^ije hadronization phase as well as the intermediate phase of constituent quarks, is proposed. For this 
•Sbdel the self-consistent solutions, which take into account the fermion vacuum polarization effects, are 
^fc^nd in (1+1) D. The renormalized total energy of the bag is studied as a function of its geometry and 
>ological (baryon) number. It is shown that in the case of non-zero topological charge there exists a set 
> configurations being the local minima of the total energy of the bag and containing all the three phases, 
■"!5^ile in the non-topological case the minimum of the total energy of the bag corresponds to vanishing size 
the phase of asymptotic freedom. 
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1 Introduction 



At present the most perspective approach to the description of the low-energy structure of baryons is 
reahzed by the hybrid chiral models (HCM) of quark bags [1-3]. The inital formulation of HCM is given 
in terms of the little bag with massless quarks and ghions confined in the phase of asymptotic freedom, 
surrounded by the colorless purely mcsonic phase, described by some nonlinear theory like the Skyrme 
model [3-5], while boundary conditions between the phases represent the "chiral confinement'' [6]. In (1+1) 
D these conditions coincide with the exact bosonisation relations, what is the origin of assumption, that 
the fermion theory inside the bag and the meson one outside are actually equivalent [7]. As a consequence, 
all the physical properties of the bag should be independent of the choice of the boundary surface, what 
is the essence of the Cheshire Cat Principle (CCP) [7,8]. However, both the exact bosonisation and so 
the CCP can be rigorously proved only in (1-1-1) D, whereas in (3-1-1) D even the formulation of such a 
problem is ambiguous. As a result, in the (3+1) dimensional HCM based on the CCP there exists a rather 
small set of observables (e.g. the topological charge), which really do not depend on the bag radius [3,8]. 

At the same time, the phenomenology of strong interactions predicts unambiguously the existence of 
the characteristic confinement scale about 0.5 fm, and so in the realistic (3+1) D models the CCP should 
be strongly violated regardless on the proof of bosonisation. Moreover, in such 2-phase HCM there is no 
place for massive constituent quarks, whose concept has been shown to be very efficient in the hadronic 
spectroscopy [9]. From this point of view the most attractive situation is the one, where the initially free, 
almost massless current quarks transmute firstly into "dressed" due to interaction massive constituent 
quarks carrying the same quantum numbers of color, flavor and spin, and only afterwards there emerges 
the purely mesonic colorless phase. The first step to such version of the bag is given by the 3-phase chiral 
model, where the additional intermediate phase of interacting quarks and mesons with non-zero radial size 
is introduced [10,11]. This model allows to take self-consistently into account: i) the phase of asymptotic 
freedom with free massless quarks: ii) the phase of constituent quarks, which acquire an effective mass due 
to the chirally invariant interaction with the meson fields in the intermediate region of finite size; iii) the 
hadronization phase, where the quark degrees of freedom are completely suppressed, while the nonlinear 
dynamics of meson fields leads to the appearance of the c-number boson condensate in the form of a 
classical soliton solution, which provides the topological nature of the model as well as the corresponding 
quantum numbers. 

It is worth-while to note, that the direct quark-meson interaction is also considered in a number 
of other approaches to the description of low-energy hadron structure, in particular, in the cloudy bag 
models [12-14], as well as in various versions of the chiral quark-soliton models [15-19]. However, the 
role of this interaction in each of these models is substantionally different. In the cloudy bag models 
such Trgg-coupling is considered only perturbatively, while in quark-soliton models it is the main nonlinear 
mechanism of dynamical generation of the quark bag in the whole space. In the case under consideration 
an intermediate variant is realized, where the contribution of the direct chiral quark- meson coupling to 
the properties of the system is nonlinear, but the confinement of quarks is provided by some additional 
procedure. Such an approach allows to realize the nonlinear mechanism of dynamical mass generation in 
the intermediate region, but unlike the quark-soliton models doesn't have problems related to the absence 
of the total confinement. 

In the present paper a toy (1+1) D model of this kind is considered, in which in the intermediate region 
the one-flavor fermion field is coupled in a chirally invariant way to the real scalar field, which possesses a 
nonlinear soliton solution in the exterior region. For this model the self-consistent solutions with different 
values of topological charge are found by taking into account the effects of fermion vacuum polarization. 
Within these solutions the renormalized total energy of the bag is studied as the function of its geometry 
and the topological charge of the solution. It is shown that for non-zero topological charge there exists 
a set of configurations being the local minima of the total energy of the bag and containing all the three 
phases, while in the nontopological case the minimum of the bag's energy corresponds to vanishing size of 
the phase of asymptotic freedom. 
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2 Lagrangian and equations of motion 



The division of space into separate phases is performed by means of the system of subsidiary fields 6{x). 
To describe the underlying machinery, let us consider the Lagrangian of the form 

Co = k {d,^4>f - ov{(t>) + \ {d^ef - glw{e), (2.1) 

where the coupling constant of sclf-intcraction of the field 9 is assumed to be large enough to neglect 
the matter fields (j) in the dynamics of 9 to the leading order, and thereafter to use 9 as background fields 
for the dynamics of 0's [10,20]. One can obviously insert in (2.1) as many fields 9{x) as needed with 
the appropriate self-interaction which will determine (almost) rectangular division of space into regions, 
corresponding to different phases, while the Lorcntz-covariance will be broken only spontaneously, on the 
level of solutions of equations of motion. Therefore in order to restore the covariance one can use freely 
the framework of covariant group variables [21]. Assuming further that subsidiary fields 9{x) have already 
formed the required bag configuration, let us start with the following Lagrangian: 

C = ^id^ + \ {d^^f - \ ml<p% - y e^^^^'^v] _ On - [V^, e'^^^'^V'] _ + V{<p)^ 9ni , (2.2) 

with 9i = 9(\x\ < xi) , 9ii = 9{xi < \x\ < X2) , 9iii = 9{\x\ > X2) being the step functions, selecting 
the corresponding regions of the bag. The commutator over fermion fields in terms, corresponding to the 
chiral fermion-boson coupling, ensures the charge conjugation symmetry of the model. Note that in this 
model the vacuum pressure term appears to be redundant, since due to the existence of the intermediate 
phase the Dirac's sea polarization behaves very specifically and itself yields the required "inward pressure" . 
Moreover, there is no special need in "valence" quarks, since it is the boson condensate in the form of the 
topological soliton which accounts now for the topological quantum numbers of the bag. 

So in initial stage we have the theory of two fields, the fermion and the boson ip. In the region / the 
bosons have the mass mo, in the region // the field (p is massless and interacts with fermions in a chirally 
invariant way, what results in emergence of the effective fermion mass M. In the region III the effective 
fermion mass becomes equal to Mq and the self-interaction of switches on, providing the appearance of a 
soliton solution for the boson field in this region. To form the bag we assume the mass Mq to be very large, 
what leads to the dynamical suppression of the fermion field in the exterior region. Proceeding further this 
way, we take the mass mo infinitely large, then for convergence of the bag's energy the boson field should 
vanish in the region /. We'll assume also that the self-interaction potential V{(p) is a an even function. 
Then the solution of equations of motion for the boson field could be either odd (the topological charge is 
by convention equal to one), or even (the topological charge is zero) function. 

Let us consider now the behavior of fields in more detail. According to the general approach accepted 
in hybrid models, wc consider the boson field in the mean-field approximation, i.e. it is assumed to be a 
c-number field. Neglecting temporarily the explicitly Lorentz-covariant description, we will consider the 
center of mass system of the bag, where p{x) should be a stationary classical field being a background for 
evolution of fermions. In the region I (p{x) = 0, while in the region III it decouples from fermions due to 
the infinite effective mass of the latters and is formed uniquely by the self- interaction V{(f). 

The equations of motion read: 
in the region / 

id^ = , (2.3a) 



in the region II 



and in the region /// 



(p = , (2.36) 

[id - Me^s^sv^ V = , (2.4a) 

[id - Moe'^^svj V' = , (2.5a) 

+ V'{^) = , (2.56) 
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where ( ) in Eq (2.4b) stands for the expectation vahic with respect to the fermionic state of the bag. To 
simpHfy calculations, we put further g = 1, because the dependence on it can be easily restored by means 
of the substitution (p ^/g- Then the spectral problem for fermionic wavefunctions V'u) with definite 
energy to reads 

c^V-o; = -ia^l + /3e'^^^ [MOji + Mo9ni] V-o, ■ (2.6) 

Upon taking Mq — oo, we get that — in the region III in such a way, that the term M^tp^ in eq.(2.6) 
vanishes, and the following boundary conditions at the points ibx2 appear 

±ij^^U^X2) + e^'^s'^(±=^2)^^(±X2) = . (2.7) 

Supplied by the condition of continuity for ip{x) on boundaries between the regions I and //, the eqs.(2.7) 
provide the correct formulation of the spectral problem for fermions, confined in the bag. Note, that the 
boundary conditions (2.7) are actually the standard chiral boundary conditions for the hybrid models [1-8]. 
However, they arise now as a direct consequence of an infinite mass of fermions in the region ///, rather 
than from a local surface action, what is not completely correct [10,22]. In the region I the Eq.(2.6) is the 
equation for free massless fermions 

ui/jj = -iaijjj , (2.8) 

while in the intermediate region II one has 

u^jj = -iai/j'jj + pMe''^''^ij:n . (2.9) 

Conditions of wavefunction's continuity on the boundary between I and II read 

i;i{±xi) = i^iii^xi) , (2.10) 

while at points |a;| = X2 the wave functions satisfy the boundary conditions (2.7). Meanwhile the field tp in 
the Eq.(2.9) is not arbitrary but has to be determined self-consistently from Eq.(2.4) with corresponding 
continuity conditions at points \x\ = x\^2- 



3 Solutions with non-zero topological charge 

The essential feature of this bag configuration is the fact that the coupled equations (2.4) in the closed 
intermediate region // of finite size d = X2 — xi possess simple and physically meaningful solution which 
would be unacceptable if these equations were considered in the infinite space. In order to obtain this 
solution in the most consistent way, we perform firstly in the region // the chiral Skyrme rotation 

xp^ = exp(-i75V?/2)xc^, (3.1) 

whatafter the Eq. (2.9) and the boundary conditions (2.7) transform correspondingly into 

{u-l ip')x^ = -iaxi + PMxo,, (3.2) 

±i-f^Xw{±X2) + xU±X2) = . (3.3) 
It follows from Eq.(3.2), that if we assume the linear behavior for the field (p{x) in the region II, namely 

if' = const = 2A , (3.4) 

then it becomes the equation for free massive fermions 

i^X = -iax' + PMx , (3.5) 

with eigenvalues u = u! — \ . So the fermions being massless in the region /, acquire the mass M in the 
region // due to the coupling to the field (p, whence the intermediate phase emerges describing massive 
quasifree "constituent quarks". 

At this moment, the most important feature of Eq.(3.5) is that it reveals the sign symmetry v — > —v, 
which corresponds to the unitary transformation of fermionic wavefunction 

X^X = niX > (3.6) 
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while the chiral currents 

j5 = i'4)-f^e'-'^^'^il) = ix^^ix (3-7) 
coincide for these sign-symmetric states: 

J5 = ix^liX = i5tliX = h ■ (3-8) 

However, the sign symmetry of Eq.(3.5) itself cannot provide the corresponding one for the fermion spec- 
trum, since it takes place in the region // only, while the latter has to be determined from the Dirac 
equation on the unification of the regions / -|- II. Meanwhile in the region / one has the Eq.(2.8), which 
possesses another symmetry, namely u ^ —to. That means that the sign symmetry u ^ —v of the 
fermionic spectrum could hold only for discrete values of the derivative in the region //. These values 
should be determined from the transcendent algebraic equation for fermionic energy levels, which is ob- 
tained from the straightforward solution of Eqs. (2.8) and (2.9) with account of boundary conditions (2.7) 
and the constraint (2.10), and reads 

1 _ „-2ikd M-i{u+k) , _ 2ikd M-i{v-k) 

exp^4za;xij ^ ^^.^^ uMiy-k) . _ oikdM+iiu+k) ' ^"^-^^ 
^ M+i(i/+fc) ^ ^ M+i{y-k) 

where = ^ . Analysing the Eq. (3.9), one easily finds, that the fermionic spectrum reveals the 
symmetry v •s-^ — z/, if 

4Aa;i = tts , (3.10) 

where s is integer, since for such values of ^'{x) in the region II the l.h.s. of Eq. (3.9) reduces to 
(— 1)* exp(4zi/xi). 

Assuming the condition (3.10) to hold, the following consequence of arguments becomes reasonable. In 
the r.h.s. of Eq.(2.4b), which determines (p"{x) in the region I/, we have the v.e.v. of the C-odd chiral 
current 

'i>,i^^^'^^t>^] =\ [x+,nix]_ , (3.11) 



with X being now the secondary quantized Dirac field in the chiral representation (3.1) 

x(a;,i) = E^"Xn(a;)e-*""*, (3.12) 

n 

where Xn{x) are the normalized solutions of the corresponding Dirac equation, 6„ , 6+ are fermionic creation- 
annihilation operators which obey the canonical anticommutation relations 

{K, = ^nn' , {bn, bn'}+ = • (3.13) 

The average over the given bag's state includes, by definition, the average over the filled sea of negative 
energy states < plus possible filled valence fermion states with cOn > which are dropped for the 
moment because their status is discussed specially below. Finally, 

(J5) = {J5)sea = k I E - E I X^illXn ■ (3.14) 
\a;„<0 u>n>oJ 

Let us emphasize here, that in Eq. (3.14) the division of fermions into sea and valence ones is made in 
correspondence with the sign of their eigen-frequencies LOn, which differ from sign-symmetric i/„ by the shift 
in A 

= J^n + A , (3.15) 

and so do not possess the sign symmetry lo <-> —to. However, if we suppose additionally, that Vn and A are 
such that for all n the signs of Un and a;„ coincide, i.e. after shifting by A none of u^s changes its sign, 
then the condition a;„^0 in Eq. (3.14) will be equivalent to the condition i^n^O . Hence 

{J5)sea =k I E - E I X^i^lXn = (3.16) 
\i/„<0 v„>o) 
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by virtue of relation (3.8). In turn, it means that (2.4b) in the region // reduces to = 0, what is in 
complete agreement with our initial assumption that (p'{x) = const in the region //. In other words, we 
obtain the solution of the coupled Eqs. (2.4) in the region // in the form 



where A takes discrete values from (3.10) and the fermion energy spectrum is determined from the relation 
(3.15), while I'n is defined from Eq. (3.9) after replacing the l.h.s. to (—1)^ cxplAiuxi). 

There are the following keypoints that make this solution meaningful. The first is the finiteness of the 
intermediate region size d, because for an infinite region // the solution (3.17) should be unacceptable. 
In our case, however, the size of the intermediate region is always finite by construction, while the boson 
field (p{x) acquires the solitonic behavior in the region III due to self-interaction V{(p). Here the following 
circumstance manifests again: in (1+1)D the chiral coupling ilje'^^-''-^ip itself cannot cause the solitonic 
behavior of the scalar field by virtue of the effects of fermion-vacuum polarization only, i.e. without 
additional self- interaction of bosons [23]. The second point is the discreteness and the u —u symmetry 
of the fermionic spectrum, what leads in turn to a reasonable method of calculation for the average of the 
chiral current J5 over filled Dirac's sea (3.16), as well as for other C-odd observables like the total fermion 
number. After all, in the case we consider the boson field is continuous everywhere and so is topologically 
equivalent to that odd soliton, which would take place in absence of fermions due to the self-interaction 
V{(p) only. That's why the topological number of the boson field doesn't depend on the existence and 
sizes of the spatial regions containing fermions (the regions / and //). On the other hand, the baryon 
number of the hybrid bag is, by definition, the sum of the topological charge of the boson soliton and the 
fermion number of the bag interior. In our case the latter is zero, hence the baryon number of the bag is 
determined by the topological charge of the boson field only and doesn't depend on the sizes of the regions 
/ and // containing fermions, what is in agreement with the general ideology of hybrid models. More 
detailed discussion of this solution of Eq.(2.4) and the arguments in favor of its uniqueness are given in 



It is also worth-while noticing, that although the (topological) quantum numbers of such a bag are 
determined by its solitonic component, it doesn't mean that the filled fermion levels with positive energy 
shouldn't exist at all. This could take place for small enough values of the parameter A only. If A increases, 
the negative levels a;„ = — -I- A will inevitably move into the positive part of the spectrum. The change 
of sign of each such level will decrease {Q)sea by one unit of charge, but if we fill the emerging positive 
level with the valence fermion, then the sum Qyai -j- Q sea remains unchanged. Analogously, the total axial 
current will be equal to J^ai + Jsea and won't change either, what ensures the vanishing r.h.s. of Eq. (2.4b) 
and so preserves the status of linear function (3.17) as the self-consistent solution of the field equations. 
Therefore, the existence or absence of valence fermions in such construction of the ground state of the bag 
depends actually on the relation between A and | 

i^\min and so appears to be a dynamical quantity like the 
other parameters (the size and mass) which are determined from the total energy minimization procedure. 

Another essential feature of this bag configuration is that (3.17) provides the self-consistent solution 
of Eqs. (2.4) for even values s = 2r in (3.10) only. The reason is that for odd values s = 2r + 1 the 
fermionic spectrum obtained from the solution of Eqs. (2. 7-10) under condition f' = 2A will always contain 
the nondegenerate energy level Xo(,x) with zero frequency vo = 0, while for even values s = 2r all the 
f„ 7^ 0. According to the general theory [24], such a zero mode yields fractionalization, what means, that 
its contribution to all C-odd observables is given by the operator ^ {b^bo — ^^o^o") with the eigenvalues ±1/2 
and the numeric coefficient determined by Xoi^)- Now let us note, that in Eq.(2.4b) the chiral current 
should be averaged over its eigenvector in order to provide the vanishing dispersion of the r.h.s., otherwise 
the system of Eqs. (2.4) would be ill-defined. So the operator part of the zero mode contribution to the r.h.s. 
of (2.4b) turns unavoidably into the factor ±1/2, meanwhile the wavefunction Xoix) appears to be such 
that the corresponding chiral current in the region // does not vanish (it is proportional to exp(— 2M|x|)). 
Hence the r.h.s. of Eq. (2.4b) doesn't vanish for odd values s = 2r + 1, and the function (3.17) is no longer 
the self-consistent solution of Eqs. (2.4). 

It's not difficult, however, to find the way of constructing analogous bags, where the odd values s = 2r+l 
are allowed. This method is based on a specific for such two-dimensional bag models possibility to choose 
in the model Lagrangian independently the signs of the chiral fermionic masses M, Mq, both to the right 




2A(x — xi) , xi < X < X2 , 
2A(x + xi) , —X2 < X < —xi , 



(3.17) 



Ref.[10]. 
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and to the left of the region of asymptotic freedom. More specifically, the considered configurations with 
even s appear from the Lagrangian (2.2), which is symmetric in signs of fermion masses to the left and to 
the right. Now let us consider the situation, when in the Lagrangian (2.2) all the chiral masses to the left 
of the central region change their signs. As a result, we obtain the following Lagrangian 



Mo 
2 



(3.18) 



where 9i = 6{\x\ < xi) , O^ff' = d{xi < ±x < X2) , of^] = ^(±a; > X2) , which is antisymmetric in signs of 
the fermionic masses to the left/right. Qualitively their behaviour for these two cases is shown on the Fig. 
1. 
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Fig.l. The behaviour of the chiral fermionic masses for the symmetric (up) and antisymmetric (down) 
cases. 

It's easy to see, that the change of the sign of Mq to the left leads merely to the change of the sign in the 
corresponding boundary condition. Namely, in this case we will have at the point x = —X2 

i-f^M-X2) + e^^«'^(-^2)Va;(-X2) = , (3.19) 

instead of (2.7). In the intermediate region after the Skyrme rotation the Dirac equation (3.5) turns into 

ux = -iax' ± PMx , (3.20) 

where in the r.h.s. of (3.20) the sign it corresponds to the sign of the chiral interaction to the right/left. In 

this case again the i' ^ —u symmetry of the spectrum takes place for the transformations X ^ X = ^2X1 
which do not alter the chiral current J5 = J5, and so the solution based on the linear Ansatz (3.4) and 
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CO = v + X remains valid, while the spectral equation for fermionic eigenfrequencies appears firstly in the 
form 

1 _ „-2ikd M-i{v+k) I _ 2ikd M-i(u-k) 
(A- \ M-i{u-k) ^ ^ M-i{u+k) 

exp(4za;xi) = - ,,,, m+.(.-.) ^ ^.... M+.(.+.) ' (3-21) 

M+i(u+k) ^ ^ M+i{v-k) 

and differs from the Eq. (3.9) by the sign in the r.h.s. The origin of this additional sign is the antisymmetry 
of the Lagrangian (3.18) in signs of the fermionic mass terms. The sign symmetry of the spectrum v <-> —v 
is still ensured by the condition (3.10), but now the zero mode emerges for even values of s, hence the 
allowed for the validity of the linear solution (3.17) values of s are the odd ones, for which exp(4iti;a;i) = 
— exp(4ii/a;i). As a result, both in even and odd cases the final spectral equation for fermions is the same 
one, namely 

1 _ „-2ikd M-i{u+k) , _ 2ikd M-i{u-k) 

exp^4zi/xi; - ^ 2ikdM+ii!^ . _ oikd M+iju+k) ■ y-^-^^) 

^ ^ M+i{iy+k) ^ ^ M+i{u-k) 

It might seem, that the odd case is physically unacceptable, since the Lagrangian (3.18) is not invariant 
with respect to spatial reflection. However, it is caused by that circumstance, that the considered models 
contain the minimal set of terms just to describe a single isolated bag. It's easy to see that the model, 
containing two odd bags, which are the mirror copies of each other, should be P-invariant. Therefore, both 
even and odd bags possess actually the same physical status, and together cover all possible values of s in 
Eq.(3.10) with the same spectral eq.(3.22) for the eigenfrequencies Vn- In other words, there exists a specific 
discrete symmetry in such systems, since the most essential properties of the bag, e.g. the confinement of 
fermions, properties of the intermediate region, the spectrum of i^^'s, etc., are invariant with respect to 
this change of signs in the Lagrangian, which predicts merely the allowed values of s in Eq.(3.10). 

Note also, that the Fig.l does not present all the possible types of the bags which could be obtained 
combining the signs of the chiral masses in all possible ways. In particular, there exist such combinations 
for which the fermions should possess an energy spectrum, which is different from that of Eq.(3.9), and 
the corresponding bags — reveal a set of different physical properties. The study of this question will be 
presented separately [25]. 



4 The total energy of the bag for the non-zero topological charge 

Without loss of generality we'll consider the bag configurations with even values s = 2r, which take place 
in the P-invariant model (2.2). The resulting configuration of the boson field reveals the structure shown 
on Fig.2. 




Fig.2. The configuration of the boson field for an isolated bag with non-zero topological charge. 
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In the region / = 0, in the region II the boson field is the hnear function (3.17), which after restoring 
the (/-dependence is sewn together with the solution of Eq.(2.5b), which is responsible for the dynamics 
of bosons in the bag's exterior. To avoid the dependence on the structure of self-interaction V(ip), we'll 
suppose, that in the region /// the asymptotic expansion of the soliton solution of Eq.(2.5b) for large |x| 
can be used, namely 

<Psolix) = ^{l-Ae-"''^) , x>X2, (4.1) 

with m being the (nonzero) meson mass in the bag's exterior, while for x < —X2 <Psol{x) is determined by 
oddness. The chiral symmetry in the exterior region of the bag is obviously lost. It should be noticed, 
however, that it occures due to specific features of (1+1) D field theories, which make the presence of the 
meson mass the necessary condition for the required soliton profile to be formed. 

The factor -k/q means that we deal actually with the phase soliton with the total amplitude being 
multiple of 2TT/g, since it is a period of the initial chiral interaction ijj exp{ij5g{p)tp. The constant A is 
determined from the continuity conditions for boson field at points x = ±X2, what gives 

xi = ^{x2 + l/m) , d = T'/"' . (4.2) 
r + 1 r + 1 

The condition d > yields then an additional restriction for the size of the confinement region , i.e. for 
the external size X2 of the bag 

mx2 > r , 

which shows, that r could be naturally interpreted as the index enumerating the excited states of the bag, 
whose sizes increase with r. 

Proceeding further, we find the relation between the parameter A and the bag's size: 

2A = Trm , (4.3) 

171X2 + 1 

whence the total energy of the bosonic soliton can be represented as 

TT^ r + 1 

E<P = m-^—. (4.4) 

g"^ mx2 + 1 

The total energy of the bag is the sum of E^p and of the fermionic contribution 

Ebag = Etp + E^. (4.5) 

As it follows from (4.4), the boson field energy decreases smoothly for increasing X2 without producing any 
vacuum pressure, despite of the fact that in the region II the gradient of ip gives rise to the constant positive 
contribution to the energy density ^ (p'"^ = 2A^, what could be identified with the vacuum pressure B in 
the standard HCM. Actually, it is an artifact of one spatial dimension in our problem: when the bag's size 
increases, the gradient of f in the region II decreases equally in any number of spatial dimensions, while 
the volume of the region // in one space dimension increases only linearly and so cannot compensate the 
decrease of A, what would take place in 2- and 3- (space) D. Thus, in (H-l)D all the non-trivial dependence 
of the total bag energy E^ag on the model parameters could originate from the fermionic contribution E^ 
only, which is the sum of the filled Dirac's sea of negative energy states and positive energy valence fermions 

E^ = Eyal + Esea ■ (4.6) 

For the ground state of the bag described above, the sum (4.6) can be reduced to a single universal 
expression by taking into account, that the charge conjugation symmetry dictates the following definition 
of the Dirac's sea energy [23, 26] 

Esea = \ - I ^ ^^n ■ (4.7) 

If the transformation from Un to Vn is sign preserving for all n and so there are no valence fermions in the 
ground state of the bag (to provide the vanishing v.e.v for the charge and chiral current), one finds from 
(4.7): 

E^ = Esea = k Yl (-^" + ^) - 5 (i^n + A) = - X l/„ . (4.8) 

t'n>0 Un>0 Un>0 
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Note, that in (4.8) the inequahty is strict because there are no levels with z/ = 0. 

If the parameter A appears to be large enough, so that the initially negative level Un = — |fn| + A 
changes its sign, it turns into the filled valence state. The latter is again necessary to provide the vanishing 
v.e.v for the charge and axial current. In this case it is convenient to calculate E^f, in two steps. First, we 
consider the contribution from all states with li^^l > A to Egea, which in analogy to (4.8) reads 

Kea = - E • (4-9) 

To this expression the energy of emerging valence fermions E^ai = — Wn\ + ^ and the contribution of the 
positive levels with a;„ = itlt'^l + A to the Dirac's sea energy should be added, what yields 

= -l^nl + A - i + A) + {\l^n\ + A)] + £^La = ~ J2 , (4-10) 

i/„>0 

i.e. the same expression (4.8) as we have got for the energy of fermions without filled valence states. 

For what follows it is convenient to introduce a set of new parameters, in terms of which the total 
energy of the bag will be expressed in the most appropriate form. First, we introduce the dimensionless 
quantities 

a = 2Mxi , 13 = 2Md, p = 2Mx2 , (4.11) 

and consider in more detail the Eq.(3.9), which determines the energy levels f„. This equation has two 
branches of roots. The first one corresponds to real k and in terms of parameters a and /3 can be transformed 
into the following form 



tan 



/ r— — ^\ X a;cos/3x + sin/3x -,r>\ 
{ayl + x^] = I = , (4-12) 

^ ' V,T^ + 1 1 — COS px + x sin px 



where the unknown quantity is the dimensionless x defined through k = Mx, so that u = M\/l + x'^. These 
real roots Xn belong to the half axis < Xn < oo, since the fermionic wavefunctions are in fact the standing 
waves in a finite spatial box with degeneracy in the sign of k, while the corresponding frequencies Vn lie in the 
interval M < Un < co. The second branch corresponds to imaginary k = iMx, u = My/l — x"^ , < a; < 1 
and can be derived from (4.12) by means of the analytical continuation 



tan 



/ I X X cosh f3x + smh I3x , , . 

ay I — x^] = , — — — . , ^ . (4.13) 

V ^ V 1 — a; cosh px + X sinh px — 1 



For this branch < i^n < M. 

Therefore, Un and so E^ appear to be functions of two dimensionless parameters a and /3, which are 
not independent, however, but give in the sum the dimensionless total size of the confinement domain p: 

a + (5 = p . (4.14) 

Proceeding further, it is convenient to extract the mass of the "constituent quark" M from the sea energy 
and fermionic frequencies as a dimensional factor: 

En = vjM = ^Jl + xl , (4.15) 

hence E^ = —M ^n- Upon introducing the dimensionless ratio of the two mass parameters of the model 

H = m/2M , (4.16) 

the dimensionless energy of fermions = E^/M and analogously the dimensionless total energy Shag = 
Efjag/M, for the latter one finds: 

5 MP + 1 

where the dimensionless parameters a, (3 are defined directly from p, and p 

« = ;:^(P + V/^), /3 = ^^- (4-18) 
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So the total energy of the bag depends ultimately on two dimcnsionlcss parameters — fi and p, where the 
parameter /7, is fixed by the ratio of the masses m and M, while the optimal value of the bag's size should 
be found from the condition of minimum of the total energy Sbagip) for given /x. 

To study the behaviour of £hag{p)^ first of all we have to renormalize the fermion sea energy which 
obviously diverges in the upper limit. Let us start with the asymptotics of roots of Eq.(4.12) in the UV 
domain, when x„ S> 1. Representing the Eq.(4.12) as 

sinaVl + = ^ (VT+~x^ + x) sin ^a-\/l + x'^ + Px + + ^ (VT+~x^ — a;) sin (^a^/l + x'^ — Px — , 

(4.19) 

where 6 = arctan x , one finds that 

p -it/2 + im 

In the expression (4.20) the first term yields the quadratic and linear divergences in and the second 

one produces the logarithmic one, while the term with the sine doesn't yield any divergence at all. To 
compensate the contribution of the first term, the energy of the sea of free fermions contained in the same 
"volume " p should be subtracted, while the logarithmic divergence proportional to P/2, is compensated 
by the relevant one- loop counterterm of the boson self-energy [10]. The remaining logarithmic divergence 
corresponding to the term l/(7r/2 + vrn) doesn't depend on the bag parameters, and originates from the 
fermion confinement inside the bag, rather than by some local interaction. Actually, it is the (infinite) 
energy of interaction between fermions and the confining potential (bag boundaries). The appearance of 
such diverging surface energy in £^ is a specific feature of fermion vacuum polarization in all the bag 
models [3,8,27-32]. 

In the considered 3-phase bag model this effect acquires some additional features. First, it takes place 
for nonzero size c? 7^ of the intermediate phase only, while the emerging surface energy is negative and 
diverges as (— 1/{tt/2 + im)). More concretely, if a — > p then (— 1)"'+-'^ sin [(7r/2 + Tm)a/ p] — —1, hence 
there remains only the logarythmic term /3/(vr + 27rn) in the asymptotics (4.20). Therefore in this limit 
£^ becomes finite just after subtraction of the energy of perturbative vacuum and adding the one-loop 
counterterm. On the other hand, the limit a — > p is equivalent to /3/a — > 0, and so the infinite interaction 
energy between fermions and bag boundaries takes place for d 7^ and the finite size of the central region 
(of the phase of asymptotic freedom) of the bag only. 

So the considered 3-phase bag model doesn't actually reveal the ability to the smooth transition into 
a 2-phase configuration for d ^ 0, although such an opportunity exists formally on the level of the 
initial Lagrangian (2.2). In fact, in the case of the two-phase bag (d = 0) the exact values of are 

= (7r/2 -|- im)/ p, and so £^ becomes finite after a single subtraction of the energy of the perturbative 
vacuum. Therefore the transition between 2- and 3-phase bag configurations requires an infinite amount 
of energy, what is a specific feature of such many-phase systems. Note also, that the ultimate role of the 
intermediate phase is the dynamical generation of the fermion mass, whereas in the case of the two-phase 
bag {d = 0) the massless fermions are reflected directly from the bag boundaries. From this point of view, 
there is an intimate connection between the infinite surface energy of the bag and the circumstance, that 
for d ^ the boundaries of the bag reflect massive fermions. 

Within the considered 3-phase bag models we have an opportunity to demonstrate this effect in a even 
more apparent way. For these purposes let us consider a P-invariant model describing the H-l-dimensional 
analog of a "dibaryon", i.e. the configuration with the topological charge 2. Such an object consists 
of two identical topological bags of the type described above, which are placed so close to each other 
that their neighbouring intermediate regions overlap. Upon dropping the ^-dependence for simplicity, the 
corresponding Lagrangian reads 

{ei + em) - 

+ Viif)^ div , (4.21) 

where 9i = d{\x\ < xq) ,9)^' = 9{xq < ±x < xi) ,9ni = 9{xi < \x\ < X2) ,9iv = 9{\x\ > X2). The 
behaviour of the chiral fermionic masses for this case is shown on the Fig. 3. 
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Fig. 3. The behaviour of the chiral fermionic masses for the dibaryon configuration. 

Again, we seek for the self-consistent solution of the model (4.21) corresponding to such "dibaryon" con- 
figuration, assuming the linear behavior (3.4) for the boson field in the intermediate regions and taking 

account of the sign-symmetry u —v as well as of the conservation of the chiral current = for the 
transformations X = (^2X- Omitting some straightforward, but lengthy calculations, let us present the 
main results. 

The profile of the boson field corresponding to the dibaryon configuration is shown on the Fig.4. 




Fig.4. The boson field profile for the "dibaryon". 

For the intermediate regions of this dibaryon one obtains ip' = const = 2A , where A satisfies the condition 

2Aa = TTS , a = xi — xq . (4.22) 

The latter condition is quite analogous to Eq.(3.10) for a single isolated bag, because the parameter 2xi in 
(3.10), as well as a in (4.22), has the sense of the size of the region of asymptotic freedom of a single bag. 
Note also, that on the contrary to the single bag models considered above, in the case of the dibaryon there 
are no zero modes in the fermionic spectrum for any values of s. So there are no additional restrictions 
imposed on the integer s in Eq.(4.22). 

It should be clear, that the H-l-dimensional model (4.21) cannot be considered as a realistic model of 
the dibaryon to any extent. However, being simple and non-trivial simultaneously, it turns out to be a 
very ftiiitful illustration for the study of the origin of additional logarithmically divergent terms l/(7r/2 + 
7rn) in the UV-asymptotics of the fermionic spectrum in such 3-phase bag models. The latter is again 
obtained from the corresponding transcendent equation for fermion levels, written in the trigonometric 
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form analogous to that of Eq.(4.12): 

sin ^2a\/r+~a^j ^x\/T+~x^cos ((/? + 7)0:; + S) — xcosjx^ + 

+ cos (^2aV 1 + a;2 j (^x^ sin ((/? + + (5) — + sin 7X + sin 7a: sm{(3x + 7)^ + 

+ (Vl + x^- cos{l3x + 7)) sin7x = , (4.23) 

where a = Ma , /3 = 2Md , d = X2 — xi , 7 = 2Mxo , 5 = arctanx. The parameter d = X2 — xi is the 
size of the external intermediate regions for each of the single bags forming the dibaryon, while 2xo is the 
size of their common internal intermediate region, i.e. the region of their mutual interaction. In this case 
the UV-asymptotics of £„'s has the following form 

P 

^ i-ir+^ (sin [(7r/2 + 7rn)(2a + j)/p] - sin [(7r/2 + 7rnh/p]) + ! + (/? + ^)/2 ^ ^ ^ 

7r/2 + 7rn ' 

where 

p = 2a + P + j = M{2a + 2d + 2xo) = 2Mx2 (4.25) 

is the total dimensionless bag's size. As in the case of a single isolated bag, the main divergent term 
in the asymptotics (4.24) corresponds to the sea energy of free fermions in the "volume" p, while the 
logarithmic term, proportional to (/3 + 7)/2, is exactly compensated by one-loop self-energy counterterm. 
The change of the coefficient in this term compared to (4.20) is caused by the fact, that in the considered 
case the interaction between fermions and boson field takes place in the region of the size 2d+ 2xo- Besides 
this, there remains again a logarithmically divergent term \/{tx/2 + 7rn), that corresponds to the (infinite) 
energy of the interaction between fermions and the confining potential (bag boundaries), and it follows 
from Eq.(4.24), that on the level of divergent terms the surface energy of the dibaryon coincides exactly 
with that of a single isolated bag. So we are led unambigously to the conclusion, that it is indeed the effect 
of confinement of fermions in a simply connected region, that yields the term l/(7r/2 -|- im) in Eqs.(4.20) 
and (4.24), since in the dibaryon configuration the number of boundary points is just the same as in the 
case of one isolated bag. Note also, that the direct consequence of this statement is that in (1-f-l) D the 
dibaryon configuration can't be obtained as a result of continuous fusion of two isolated bags, since when 
they are separated enough from each other, the sum of their surface energies is twice larger than that of the 
dibaryon. In other words, in (1-|-1)D the reconstruction of the bag's surface in the fusion-fission processes 
requires an infinite amount of energy. 

After all, it follows from (4.24), that for /3 — > 0, i.e. for vanishing external intermediate regions of the 
dibaryon one gets (— l)"''"-'^ sin [(7r/2 -|- 7rn)(2Q! -|- 7)/p] — > —1, what compensates the term /2 + -Kn) , 
and the infinite interaction energy between fermions and bag boundaries disappears too. So we are again 
led to the conclusion, made for a single bag by analysis of the asymptotics (4.20), that the infinite surface 
energy appears only then fermions pass through the intermediate phase just before reflection from the bag 
boundaries. 

As a result, for a 3-phase bag with ti 7^ the extraction of the finite part from £.^, consists actually 
of two separate procedures. The first one is the standard renormalization onto perturbative vacuum 
with account of the one-loop counterterm, caused by virtual fermion pairs [10]. The second one is the 
compensation of the surface energy by means of an appropriate subtraction, and both procedures suffer 
from an ambiguity in the choice of subtraction point. In the "classical" renormalization scheme, the 
uncertainty in the choice of subtraction point is cancelled by fixing the physical values for a corresponding 
number of parameters. For obvious reasons, we avoid doing that in our "toy" (1-|-1)D model, but instead 
consider the most straightforward approach to the compensation of divergences in the sum (4.8), which 
preserves the continuous dependence of the result of substraction on the model parameters. The essence 
of this approach is that we subtract from ^„ e„ another sum with the same summation index n, whose 
common term coincides exactly with the divergent part of asymptotics (4.20). The result is the finite 
quantity 

'71/2 + 7171 1+13/2 



E 



+ 



7r/2 + 7rn 



(4.26) 
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This method requires no counterterms because all the divergences are already cancelled by the subtracted 
sum. Of course, to some extent the physical meaning of such procedure is lost. It should be emphasized 
however, that it is only the (1+1)D case when the theory with coupling Cj = G'ip{a + i757r)V' is (su- 
per) renormalizable and any counterterm has explicit physical meaning. For higher space dimensions this is 
already not true and so the procedure of compensation of divergences in the energy based on (4.26) should 
not be considered as having no motivation. For more detailed discussion on the extraction of the finite 
part from the divergent Dirac's sea energy in 3+1 D HCM see Refs. [30-32]. 
Proceeding further, let us turn to the study of the total bag energy 

as a function of the parameters p. The analysis of the contribution of convergent logarithmic part from 
the sine-term in the asymptotic expression (4.20) to £^ yields the first feature of £bag- Let us transform 
this contribution to the form 

{£i.)^Ja,(i) - - ^^172 (4-28) 

and then use the well-known relation 

g(_,^n Sm[.(n+l/2)] ^ j^^^^^^^^ ^ ^^^^ ^ ^ _ 

n=0 n-t-i/z 

It is easy to see, that the sums (4.28) and (4.29) possess the similar common term, while the sum (4.28) 
diverges as (— ln(7r — z)) when z — tt. Hence, for ira/p — tt , what means /? — 0, the sum (4.29) will show 
the similar behavior, namely: 

(S^)^ (a,/3)^--ln/3, /3 ^ 0. (4.30) 

V ^ Hog TT 

Therefore, both the regularized fermion energy (4.26) and the total bag's energy reveal the logarithmic 
singularity for /3 — 0, i.e. for p — >^ r//x, what is in complete agreement with the qualitative analysis of 
the vacuum polarization effects in the 3-phase bag, performed above. In particular, for ^ both the 
intermediate phase and the infinite energy of coupling between fermions and bag boundaries disappear, 
what in the present case shows up in the unlimited increase of fermion energy, since after the renormalization 
of £^ by means of subtraction (4.26) the divergent part of the surface bag energy is included in the onset 
of energy. 

^hag will also grow for p — > 00. In this case a/p ^ r/{r + 1), so the logarithmic term (4.28) remains 
finite, what implies, that we have to deal now with the whole sum (4.26). However, the leading order 
behaviour of £^ could be evaluated from (4.26) quite effectively by using the fact that for p — *^ 00 the 
fermionic spectrum becomes quasicontinous, what allows to transform the sum over Xn into integral over 
dx. In particular, the analysis of distribution of the roots of Eq. (4.12) shows, that for this limit J2n^n is 
approximated by the following (divergent) integral: 



En ^ — [dx VT+~X^ 



I 3,2 sin ( aV 1 + ) cos ( aV 1 + x 



13 + K + a 



l + x"^ 2;2 _|_ gjn2 (ay/T+^^ VTTx^ (x^ -|- sin^ (aVT+^^ ' 

(4.31) 

For the subtracted sum in (4.26) one can easily find 

E (^-^^ + ^j^] ^ ^ fdx (x + i±^) . (4.32) 
P Tr/2 + 7mJ TT J \ px J ^ ^ 

The integrals (4.31) and (4.32) have obviously the same divergent part 

- / dx{px + l/x + l3/2x) , 

TT J 

and so their difference yields a converging integral, in agreement with the subtraction procedure. The 
leading term of the integrand in this difference, taken with the (correct) inverse sign, is P/Sirx^. Since 
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j3/p — s- l/(r + l) for p — oo, this finally leads to the emergence of the positive, proportional to p, 
contribution to if ^., and correspondingly to £bag- 

The numerical calculation confirms completely such qualitative predictions for the behavior of S^jjip) 
and £bag{p)- In the present paper such calculation has been performed for /x = 0.25, what corresponds 
approximately to the ratio m.,^/2MQ, where the constituent quark mass Mq is assumed to be equal to 300 
MeV, and for g = 1, because the energy of boson soliton doesn't have any significant influence on the main 
properties of £bag{p)- The results of £bag{p) calculation for r = 1,2,3,4,5 are depicted on the Fig. 5 and 
show, that the size and energy of the solution, determined from the minimum of Sbagip), continuously grow 
for increasing r, whereas the curvature of Sbag{p) in the minimum decreases, what provides an unambiguous 
interpretation of configurations with r > 1 as excited states of the bag. 



5 Bags with zero topological charge 

Now let us consider the bag with vanishing topological charge, what for the relevant configuration of the 
boson field should be an even one (p{x) = (p(—x). The principal difference between this case and the 
previous one is that for even cp{x) the sign symmetry lo <^ —u is an immanent feature of the spectral 
problem for fermions (2.7-10), what can be easily justified by means of the following transformation of 
fermionic wavefunctions 

'^Ljix) tp-u>{x) = ±j5ll^u>{-x) . (5.1) 

However, the corresponding axial currents are related now in the following way 

j-Ux) = -£{-x) , (5.2) 

so there is no automatic compensation between positive- and negative-frequency terms in the v.e.v. of 
J5{x). From (5.2) one can derive only the relation 

{J5{x))sea = {Jb{-x))sea , (5.3) 

what guarantees the consistence of Eq. (2.4b) with respect to parity. The direct consequence of such 
fermion properties is that the even configuration of the boson field, similar to (3.17), 

f +2\{x - Xi) , Xi < X < X2 , , . 

fy.x)-<y _2\{x + Xl), -X2<X<-Xl, 

is no longer an exact solution of Eqs. (2.4), since in this case {Jb{x))sea ^ in the region II. 

Nevertheless, the configuration (5.4) plays an important role in the study of the non-top ological case. 
First of all, for small values of g it turns out to be a rather good approximation to the precise solution. 
To argue this statement, let us note firstly, that the substitution ip = (p/g removes g from the Eq.(2.4a), 
while Eq.(2.4b) will contain g only as a coefficient in the r.h.s, namely 

^" = i9^^{[i^,l^e''''H]_) ■ (5.5) 
Proceeding further, it would be natural to assume that the potential V{(p) depends on g as 

V{^) = W{gip)/g^ , (5.6) 

where W{f) should be an even polynom to preserve the (anti)symmetry of soliton solutions. Therefore, for 
small g there emerges quite naturally an expansion in powers of g^ in the problem. Within this expansion, 
the zero-order approximation for the rescaled boson field p>{x) is the configuration (5.4) in the region 77, 
while if>{x) = {) in the region 7, and in the region 777 it is assumed to merge with the asymptotics of an 
even soliton solution 

(psoi{.x) = (l - Ae-"'^^^) , |x|>X2, (5.7) 

quite similar to the topological case. Sewing together (5.4) and (5.7) by means of the continuity conditions 
for (p and ip' yields the following relation 

Tim, 

2A = — - — - , 5.8 
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whence for the energy of the boson field one finds 



TT^m 



Note, that in the initial variables ip the dependence on g in is restored by adding the coefficient l/g"^. 

Now let us show, that the first-order 0{g'^) correction to the energy of the boson field (5.9) vanishes 
exactly for any current in the r.h.s. of Eq.(5.5), provided the asymptotics (5.7) for the boson field in 
the region /// remains valid beyond the perturbation expansion in g^, what implies, that the corrections 
caused by the r.h.s. of (5.5) disturb only the value of the parameter A. To simplify calculations, we will 
consider further only the positive half-axis. The contribution of the negative one is exactly the same. 

Prom the asymptotics (5.7) we derive 

m(p{x2) + (p'{x2) = TT , (5.10) 

while in the region /// 

^'jjj{x) = ^'(x2)e-"^(^-^2) ^ (5,n) 

which are valid beyond the y^-expansion as well. Using the virial theorem which is also relevant beyond 
the latter expansion, we obtain the following general expression for the contribution of the region III to 

/// 

Proceeding further, on account of the first-order correction from the non-zero (J^^x)) one obtains for the 
boson field in the region // 

(f{x) = 2X{x — xi) + g'^(pi{x) . (5.13) 

At the same time, it follows from the condition <f{x) = in the region / and the boundary conditions 
(5.10) that 

(pi{xi) = , m(f>i{x2) + (f>'i{x2) = . (5.14) 
Then for the boson field energy in the region II with the first 0{g^) correction one finds 

E^ii = l jdx (p'"^ = 2\(\d + g'^(pi{x2)) . (5.15) 

On the other hand, it follows in the same approximation from (5.12) and (5.13), that 

^^m = ^(A + 5Vi(^2)). (5.16) 

Returning to Eq.(5.14), one finds, that in the sum -Ej^jj + E^jjj the contribution of (pi vanishes. In 
other words, within the 5^-expansion the corrections to the leading approximation (5.9) in Eif,, caused by 
nonvanishing {J^(x)), begin from the second order 0{g'^) only. 

Next, let us note, that for fermions the y^-expansion starts from the order 0{g^). In this approximation 
the spectral problem (2.7-10) leads to the following equation for the fermionic spectrum 



exp^Aiuxi) = 



e 



-2ik+d M-i{u++k+) 1 _ 2ik-d M-i{u--k-) 

M-i(v+-k+) ^ ^ M-i{v^+k^) 

1 _ p-2ik+dM±ii!^+i:h+l 1 _ r'2ik-d ^+ii''-+k-) 
^ ^ M+i(u++k+) ^ ^ M+i(y_-k-) 



(5.17) 



where v± = uj \, v\ = k'j^ + M^. The total energy of the bag is still given by the sum (4.5), where the 
fermion energy has the following form 

u}„<0 

In (5.18), as well as in (4.8), the inequality < is strict, because for the configuration (5.4) there are 
no levels with a;„ = for any values of Xi,X2- 
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Finally, after restoring dependence on in E^p we obtain the following expression for the total energy 
of the bag 

^^"^ = 7 dri + + ^^^'^ ' ^^-^^^ 

where the first two leading terms in Ei,ag — the bosonic 0{l/g'^) and fermionic 0{g^) — are determined 
by the zero-order approximation for the boson field (5.4.7) only, while the corrections start with 0{g^) 
terms, simultaneously in the bosonic and fermionic parts of the total energy. 

The numerical calculation completely confirms such behaviour of the system for small g. On the Fig. 6 
a typical profile of the numerical solution for the boson field is shown, obtained by means of minimization 
procedure of the total energy functional for the values of g from the interval 0.01 0.5, which shows clearly 
that for such values of g the exact solution in the intermediate region is almost indistinguishable from the 
linear function (5.4). 




Fig. 6. The profile of the numerical solution for boson field, obtained by the minimization procedure of 
the bag's total energy functional for 5 = 0.1 on a spatial lattice with ~ 10^ nodes. The self-interaction of 

bosons in the region /// is taken in the form V{ip) = (vr^ — {g^pf') /2(7^. The dashed line, which merges 
almost with the solid one, corresponds to the calculation with the doubled number of nodes on the lattice, 
what indicates the accordance of the numerical solution with the continuous limit. 



Moreover, the respective simplicity of Eq. (5.17) makes it possible to analyse the fermionic spectrum in 
a semi-analytical way, what in turn allows to use the configuration (5.4,7) as a trial one for a qualitative 
study of the ground state properties for the non-topological bag for even larger values of g ~ 1. 

Thus, in further analysis of the main properties of the non-topological bag we will use the first two 
terms in the total energy (5.19), which can be found directly from the configuration (5.4,7). Recalculating 

to the dimensionless variables, introduced in (4.11,14-16), one obtains 

£ug{a,l3) = £^{a,l5) + ■ (5-20) 
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Note, that in (5.20) the sum of the parameters a, /? defines the dimensionless size of the bag p, but on the 
contrary to the previous case in all other respects they should be considered as independent ones. So the 
study of the bag's energy as a function of its geometry becomes a qualitively different problem of finding 
the two-dimensional surface £bag(a,P). 

The extraction of the finite part from £^{a, P) undergoes the same main stages as in the topological case, 
but reveals some peculiar features, caused by the independence of a and /?. After some transformations 
it's not difficult to obtain from (5.17) the UV-asymptotics of the energy levels in the following form 

e„(a, i3) = + + (-ir^cos(2Ad)sin[(vr/2 + .n)./p] + l + /3/2 ^ ^^^^^.^ _ 
p 7r/2 + 7rn 

It follows immediately from the structure of the logarithmic term in (5.21), that as in the topological case, 
the renormalizcd via asymptotics S.^. as well as Sbag' show up a logarithmic divergence (— ln/3/7r) for 
/? — > 0. Besides this, £^ and £bag increase for /3 ^ oo and finite a. Since p grows simultaneously with (5, 
this effect turns out to be quite similar to the increase of S^p and £},ag for p — oo in the topological case: 
in the UV-domain the difference between and v- vanishes and Eq.(5.17) turns into (3.9). Thereafter 
to analyse the rcnormalized £^ one may use the integral approximation (4.31,32), in which for /3 — > oo the 
main term in the integrand of £^ is positive and proportional to f3. 

The behaviour of £^ and Shag for a — ^ oo and finite /? requires special consideration, for in this case 
the logarithmic term with the sine in asymptotics (5.21) becomes significant again, but unlike the case 
of /3 ^ 0, there appears now an additional factor cos{2\d). Since 2Xd = 7r/i/?/(/i/3 + 1), the sign of this 
multiplier can be either positive or negative depending on the current value of /3. However, in this limit 
there becomes significant another effect, namely the proportional to a increase of the number of levels on 
those branches of the fermionic spectrum, which correspond to the imaginary values of k± in (5.17), where 
< i'± < M. Directly for the Eq.(5.17) this effect shows up in an intricate enough way due to the presence 
of separate branches for imaginary and k- and therefore can be analysed in detail only numerically, 
but its essence could be understood quite simply, if we neglect for a while the difference between k+ and 
k-. Then we are left with only one branch with < Un < M determined from Eq. (4.13), from which it is 
easily seen that for a — >■ oo the spectrum of energy levels belonging to this branch becomes quasicontinous 
with the interval between the levels of order ir/a, hence J^n^n over this branch can be approximated by 
the integral 



^n—>- / dx 

n Jo 



II branch 



ax ^x{(3 + .T^)) - smh{(3x + - x'^ 



Vl — cosh(/3a; + 7) — \/l — 



2 



(5.22) 



Therefore, for a ^ 00 the second branch contribution to X^n^n takes the form a/27r + finite terms 
depending on /3 only. Transforming further the subtracted sum to the integral one obtains 



a 



/V2 + .„ 1 + {-^l±m) , (5.23) 

V V P 7r/2 + 7rny tt J \p x J ' ^ ^ 

II branch 7r/2 

whence it follows that for a 00 the main terms in the subtracted sum should be 

\ -^Ina. (5.24) 

The leading, proportional to a terms in Eqs. (5.22,24) cancel each other, so after subtraction the contribu- 
tion of the second branch to the renormalized should be (l/7r)(l -|- /3/2)lna -t- finite terms. For the 
case of separate branches for A;+ k- the general features of their asymptotic behaviour for « — >^ 00 remain 
the same. As a result, after combining this asymptotics with the corresponding input of the logarithmic 
term in the UV-asymptotics (5.21), for the leading term in £^ for a — 00 one finds 

-(l-F/3/2 + cos2Ad) Ina , a ^ 00 , (5.25) 

TT 

which is definitely positive for all (5. So in this limit the bag's energy also grows, but now proportional to 
In a. 
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The numerical calculation confirms completely such qualitative behaviour of £^ and £bag- It has been 
performed for the same values of /v, = 0.25 and 51 = 1 as for the topological bags, shown on Fig. 5. The most 
specific feature of the non-topological case is that in accordance to (5.9) the boson field energy depends now 
on (3 only, while the renormalized energy of fermions decreases smoothly for a — >^ 0. A profile of the surface 
S^{a, P) is depicted on Fig. 7, where = for a = /3 = is an artifact of the chosen subtraction method. 
As a result, there isn't any non-trivial minimum in the total energy for the non-topological case at all, 
while the minimal energy is achieved by the configuration with vanishing size of the phase of asymptotic 
freedom and for finite non-zero /?, what is clearly seen from Figs. 8,9, on which the profiles of the 2D 
surfaces £i,ag{oi,P) are presented in different scales. So for the bags with zero topological (i.e. baryon) 
charge the considered three-phase model predicts that the main role should be played by the intermediate 
phase of constituent quarks, what is quite consistent with semi-phenomenological quark models of mesons 
[9,33]. 

6 Conclusion 

This work was aimed at the construction of a reasonable model of a hybrid chiral bag without special 
assumptions about the intimate relations between the fermionic and bosonic phases. Our results show 
that such a model can be actually formulated in a quite consistent fashion, and to certain extent could 
be more efficient way of description of low-cncrgy hadron physics compared to the standard HCM [1- 
8]. The main advantages of such an approach include, first of all, a more correct formulation of chiral 
boundary conditions for which all the components in the Lagrangian possess clear physical meaning, the 
existence of the intermediate phase, describing quasifree massive "constituent quarks", as well as a quite 
acceptable physically behaviour of the total bag's energy as a function of its size, which takes the form of an 
infinitely deep potential well with a distinct minimum in the topological case, whereas in non-topological 
case the minimal energy of the bag is achieved by the configuration, where the phase of asymptotic freedom 
disappears. 

Besides this, in this model the condition of fermions confinement, incorporated into it from the very 
beginning, shows up more explicitly. It manifests, in particular, in the fact that there is no need in the 
vacuum-pressure term B, which in the standard approach is inserted into the model by means of some 
extra assumptions, since in the considered case the Dirac's sea polarization itself produces the infinite 
increase of energy at large distances. Another essential feature is the appearance of infinite interaction 
energy between fermions and bag boundaries (confining potential) for d 0, what implies, that the size of 
the intermediate region doesn't actually vanish, although on the level of the initial Lagrangian the formal 
limit d ^ exists and yields the standard two-phase model of a hybrid bag. In other words, such a 3-phase 
model cannot be continously transformed into a 2-phase one, what is the ultimate reason of remarkably 
different features of this model compared to the standard 2-phase ones. 

It is worth-while to mention once more the question of the choice of method of calculation of the Dirac's 
sea averages for fermion bags. The method we used is based on the discreteness of the fermionic energy 
spectrum what by means of quite obvious considerations leads to very simple solution of self-consistent 
equations of the bag in the intermediate region. Let us remark, however, that despite of arguments in 
favour of such method of calculation of sea averages, we cannot completely reject alternative methods like 
the thermal regularization. The question of which one is more adequate to the physics of the problem, 
should be answered only by means of detailed study of realistic models. 

It should be also emphasized, that by constructing such a 3-phase model we have substantially used the 
condition of lorentz-covariance. The initial formulation of the model is a local field theory and regardless 
on the variety of classical solutions one needs to deal with, the covariance is broken only spontaneously, 
and so can be freely restored by means of methods of Refs.[21] based on the covariant group center-of-mass 
variables for a localized quantum-field system. However, such an explicitly covariant framework requires 
some essential changes in the calculation techniques, since the invariant dynamics of fields in the c.m.s. 
acquires a specific finite-difference form [21], and so will be considered separately. 

This work was supported in part by RFBR under Grant 00-15-96577 and by Sankt-Petersburg Con- 
currency Centre of Fundamental Sciences, Grant 00-0-6.2-22. 
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Fig.5. Sbagip) for r = 1, 2, 3, 4, 5 in the topological case. 
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Fig.7. The profile of the surface £^{a,P) for the non-topological bag. 
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Fig.8. The profile of the surface £ijag{(x,P) for the non-topological bag. 
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Fig.9. The profile of the surface £bag{o(,l3) for the non-topological bag, rescaled to observe the behavior 
for small /?. 
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